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ABSTRACT. We prove a remarkable formula of Ramanujan for the logarithmic deriv- 
ative of the gamma function, which converges more rapidly than classical expansions, 
and which is stated without proof in the notebooks [5]. The formula has a number 
of very interesting consequences which we derive, including an elegant hyperbolic 
summation, Ramanujan's formula for the Riemann zeta function evaluated at the 
odd positive integers, and new formulae for Euler's constant, 7. 



Ramanujan was evidently fond of series expansions and representations of spe- 
cial functions. The unorganized material in the second and third notebooks contain 
many interesting formulae involving the gamma function, Bessel functions, and hy- 
pergeometric functions, to name a few examples. The following remarkable formula 
for the logarithmic derivative of the gamma function has some very interesting con- 
sequences, and should be contrasted with the well-known series representation 

r' 00 (\ 1 \ 

i>{x + 1) := — (x + 1) = -7 + V — , 

1 \n x + n I 

71 = 1 x / 

[1, p.259, 6.3.16]. 
Theorem. For all x > 0, 
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The theorem is recorded as formula 2, p. 280 in Ramanujan's notebooks [5]. As 
per usual, Ramanujan gives no proof or explanation. In the sequel, we give a proof, 
and offer a plausible explanation of how Ramanujan may have arrived at formula 
(1). But first, we examine the formula more closely, with a view to gaining a better 
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understanding of the formula and its many interesting consequences, some of which 
may not be immediately apparent. 

The formula takes some time to absorb, but even an initial glance reveals that the 
convergence is more rapid than the classical representation. The appearance of the 
logarithm in a summand is very striking, if not mysterious, and at first may arouse 
some suspicion, since ijj(x + 1) is well-behaved for x > 0, whereas on the right we 
have terms which contribute singularities at the positive integers. However, a closer 
examination reveals that the singularities arising from the various terms actually 
cancel each other out. For example, from the partial fraction decomposition 

, 1 ^ 2x 

7rcot(7rx) = - + > — — 

x ^ x 2 - k 2 

k=l 

it can be seen that the pole contribution from 

7rcot(7nr) 

e 2irx _ I 

at x = m (m G Z + ) is cancelled by the mth term in the sum 

2k 

2_, ( e 2^ _ l)(k 2 -X 2 ) 

leaving only a harmless (finite) quantity. This is a routine exercise. Similarly, using 
the fact that 

| sin(7rx)| ~ ir\x — m\ {x — > m) 

i t i i ■ i • ■ i • •! t t vrlog |2sin(7rx)| 

one can show that the logarithmic singularity contributed by ^ & t 

2sinh (ttx) 

x = m (m £ Z + ) is cancelled by the mth term in the sum 

7T log \ k 4 — x 4 | 
2^ sinhVfc) 

Carrying out the details of the aforementioned limit calculations provides the fol- 
lowing formula for Euler's constant, 7. 

Corollary. Let m be a positive integer. Then 

2k 
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Note that for the inner sum in (2) we have [1, p. 259, 6.3.17] 

°° k 2 
^>n(k 2 + n 2 ) 

n=l 

We have used MAPLE to check (2) with various values of m. It is interesting to note 
that with m = 1, taking a mere five terms from each of the series indexed by k is suf- 
ficient to give 1— 7 accurate to thirteen decimal places: 1— 7 = 0.4227843350984 

Proof of Corollary. Let i^min (1). It suffices to show that 

, 7rcot(7rx) 2m \ 1 ir 

hm — — i — - + 



x—>m 



e 2nx _ 1 1 ( e 27rm _ i)( m 2 _ x 2) J 2m(e 2 ™ - 1) 2sinh 2 (7rm) 
and 

/vrlog |2sin(7rx)| tt log |m 4 - x 4 | \ _ 7r (log 7r - log(2m 3 )) 
x^m y 2sinh 2 (7rx) 2 sinh 2 (7rm) / 2sinh 2 (7rm) 

We have 

7r cos(7r(x — m)) 7r(l + 0(x — m) 2 ) 



TT COt(7Tx) 



sin(7r(x — m)) ir(x — m)(l + 0(x — m) 2 ) 
+ 0(x — m), 



x — m 

1 27re 27rm 



and 



g2nx ^ g27rm ^ ^g2-7rm ]^2 

2m 1 



(x — m) + 0(x — m) 
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(e 27rm - l)(m 2 - x 2 ) (e 27rm - l)(m - x) 1 + 
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1 / x — m 



(e 2nm - l)(m - x) V 2m 



1 h 0(x — m 



1 1 



(e 27rm - l)(m - x) 2m(e 2 ™ - 1) 
Thus, 

7rcot(7rx) 2m _ 1 2Tre 2nm 

e 2 ™ - 1 + (e 2 ™ - l)(m 2 - x 2 ) ~ 2m(e 2 ™ - 1) ~ (e 2 ™ - l) 2 + (X ~ 

which proves the former. For the latter, we have 

/log |2 sin(7rx)| log |m 4 — x 4 

L := lim 5 o 

x^m \ ginh (ttx) sinh (7rm) 

log 2 — log(4m 3 ) /log | sin(7rx)| log|m — x| 



+ lim 



sinh 2 (7rm) a— >m y sinh 2 (-7rx) sinh 2 (7rm) 
log(2m 3 ) /log 1 7r(m — x) + 0(m — x) 2 1 log |m — x| 



+ lim 



sinh 2 (7rm) x^m \ sinh 2 (7rx) sinh 2 (7rm) 

log7r — log(2m 3 ) /log |m — x| + 0(m — x) log|m — x 



+ lim 



sinh 2 (7rm) x— >m y sinh 2 (7rx) sinh^(7rm) 

log7r — log (2m - 
sinh 2 (7rm) 
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as required. □ 

Instead of letting x tend to a positive integer in (1), one can of course obtain a 
formula for 7 by directly substituting x = m — 1/2 (m G Z + ) in (1). Jonathan 
Borwein has used this approach to compute 7, choosing large m to balance the error 
terms. However, he reports that Brent's method is computationally superior. Part 
of the difficulty stems from the need for full precision calculation of the logarithms, 
and of e 7rk . If there is an efficient way of doing these calculations, it is possible 
that an algorithm based on Ramanujan's formula could beat other methods of 
computing 7. [J. Borwein, personal communication.] 

Another very interesting, and perhaps unexpected consequence of (1) is obtained 
by letting x — > +00. Since this limiting case is an essential ingredient in the proof 
of (1), we shall temporarily content ourselves with heuristic reasoning, deferring 
the rigorous details of the argument to the sequel. As x — > +00, most terms drop 
out, and since ijj(x + 1) ~ logx, what remains is 

00 , ,,4 4, 
7T . 7T log \k — X I . , 

logx ~ - logx- - 2^ . . 2 , (x^+00). (3) 

o 2 ^ sinh [irk) 

For large x and small k, log \k 4 — x 4 | ~ 4 logx, whereas when k is large, sinh 2 (7rfc) 
is very large, so that (3) becomes 

x 1 

logx~ -logx- 2?r logx V ; „ (x^+00). (4) 

3 ^ smh (irk) 

For (4) to be true, the coefficients of log x must agree i.e. we must have 

1 _1 1 (5) 

^ sinh 2 (vrA;) 6 2vr' y ' 

a remarkable result, with an interesting history in and of itself. See [2, Prop. 2.26]. 
That (5) is a special case of (1) makes Ramanujan's formula all the more intriguing. 

An additional consequence of (1) is Ramanujan's remarkable formula for the 
Riemann zeta function, evaluated at the odd positive integers. 

Corollary. Let N be a positive integer. Then 



U-2N-1 U-2N 

+ (l + (-l) N )n, 

k=i - k~i smh ( 7rfe ) 



2N((2N + 1) + -^—^ + (1 + (-in^T 



(2ir) 2N+1 s T(-l) j+1 (2i - 1) B2j B2N + 2 ~ 2 i (q) 
[ ' } [J ) (2j)\(2N + 2-2jy: W 



Here, B = 1, Bi = —1/2, B 2 = 1/6, etc. are the Bernoulli numbers, defined by 



oo 



k=l 
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Since the Bernoulli numbers are all rational, (6) shows that C(2^V + 1) is a rational 
multiple of 7r 2Ar+1 plus either one or two rapidly convergent series, the number 
depending on the parity of N. 

Proof of Corollary. The main idea stems from the observation that ip(x + 1) is 
essentially the ordinary generating function for the sequence C(2), C(3), • • • • More 
precisely, 

oo 

^( x + l) = - 7 + ^ x "(_l)»+i C (n+l), (|x|<l). (8) 

n = l 

Thus, we are tempted to expand both sides of (1) into power series in x and then 
equate coefficients of x 2N . However, since logx has no Maclaurin series, we first 
differentiate both sides of (1) with respect to x. Since the machinations involved in 
this differentiation arise in our proof of (1), we shall defer the details to the sequel 
and content ourselves for the moment with merely reporting that after simplifying 
the resulting expression, we obtain 

,/. . 7T 1 1 7r 2 csc 2 (7rx) 



3x 2x 2 2ttx 3 e 2 ™ - 1 



4kx ^ 2ttx 3 



E<iKX \ -> Ll\ X . . 

^ (e 2 ^-l)(fc 2 -x 2 ) 2 + ^ sinh 2 (7rfc)(fc 4 -x 4 ) • ( } 

We now expand all terms of (9) into power series in x. We have 

OO . 7 OO, OO oo 

EAkx x - 4x x - /x\ 2n - 2 

k — 1 k — 1 tl — 1 

00 00 £.-2n-l 

=*£"*"-' £5^, ( 10 > 

n=l k=l 



E° 2irx 3 ^ 27tx 3 ^ /x\ 4ra 

A . =1 sinh 2 (vrA;)(A: 4 -x 4 ) ~ ^ /c 4 sinh 2 (7rfe) ^ U J 

^ ^sinh 2 (7r£:) 

n=l fc=l v ' 

and from (8), 

oo 

V>'(x + 1) = ^(-x) n -V(n+ 1). (12) 

71=1 

Finally, using (7) and the derivative of the well-known power series for 7rcot(7rx), 
we have 

2 2 / \ 1 00 7->°° 7-) 

TTCSC^TTx) 1 V~Vo_„.\ti ™ V~V ^fc+l/oi. 1 \™2fc /r)_\2/c 2fc 



E( 2 ^) n 3 L E(- 1 ) fe+1 ( 2fc - 1 ^ 2fc ( 27r 



e 27rx _ ! 2vrx 3 ^ v ; n! ^ v 7 v y v y (2fc)! 

n=0 fc=0 V ; 
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oo n R R 

+ 27 rx^ [Z7TX) ^ { L) [ZJ ij (2j)! (2n + l-2j)! 

n=0 j = K J/ \ J) 



i i 



,2ti+1/ 1\7I+1 



(-l) n+1 (2n- 1) 



(2j)!(2n-2j)! 

-E>2n B]_ 



n=l 



(2n)! 1! ' 



(13) 



since S 2 n+i-2j = for j = 0, 1, . . . , n — 1. Substituting (10) through (13) into (9) 
yields 

oo oo oo , _2 n — 1 oo oo , —An 

y (-x)"- x nc(n + 1) =4 y nx 2 "- 1 y 4he — + 2^ y x 4 "- 1 y — % — 

i i fi e - 1 i fl sinh 2 (7rfe) 

n=l 71=1 fc = l ri=l fc=l v ' 



^ OO 71 



71 = 2 
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(2j)\ (2n-2j)\ 



27TX 3 
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The corollary (6) now follows on equating coefficients of x 2N 1 in (14). □ 

Remark 1. Equating coefficients of x 2N in (14) gives Euler's celebrated formula 
C(2AT) = { -l)N {2it? n ^_ Bi = ^n+i^n-i^nJ^^ which ig also vaM 

when N = 0, since ("(0) = —1/2 = B\. It should also be noted that (6) holds for 
N = if we make use of the hyperbolic summation (5) and interpret the term 
2N(,(2N + 1) as lim 2x((2x + 1) = 1. In fact, it can be shown that (6) holds for 

x— >0 

negative integers N as well. For example, the case N = —1 in (6) gives 



E 

k=i 



k 



1 1 



e 2^k_i 24 8tt' 



As a further example, let m > 1 be odd, and let N 
well-known evaluation ("(1 — 2m) = — - — , we obtain 



-m. From (6) and the 



2m 

00 1.2(71-1 

Efv 
e 27rfc _ I 4 m 

k = l 



B 



2771 



Recalling that 



OO 2,2771-1 



dx 



D 



2771 



, we derive a rather striking equality be- 



; e*"^ - 1 4m 
tween a sum and an integral in which the summand and the integrand are formally 
identical: 

00 fc2m-l 

2 P 2,k _ 1 = /„ e Lx _ 1 dX ' (° dd m > 



fc = l 



1 Jo 



00 2,2th-! 
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Remark 2. Let a, (3 > with a(3 = ir 2 , and let N be a positive integer. One can 
derive the slightly more general formula 

2^0"^ C(2iV + l) + 2^ 



e 2ak — \ I 
k = l I 

la'-T^ (-0) 1 - N Y 



^ sinh (afe) ^ sinh {(3k) 

N+1 R R 

= E(-D' +1 w - v 1 -'* ^wfi^y. <15) 

. , . , , , , ,. vr0racsc 2 (2:^/7ra) . 

in the same manner as our proof of (bj by expanding ■= into par- 

tial fractions and then comparing coefficients of like powers of z. The author is 
currently investigating the feasibility of this relatively elementary approach to de- 
riving analogous Ramanujan-like formulae for the Dirichlet L-functions via related 
partial fraction expansions. 

Other proofs of (15) have been given via Mellin transforms, the Poisson sum- 
mation formula, and elliptic modular transformations, to name some examples of 
various methods employed. See [2] and [3, p. 276] for an extensive bibliography. 
Since Ramanujan was quite fond of partial fraction expansions, it is possible that 
our derivation is close to Ramanujan's. 

We have seen (2) that specializing x to positive integer values in (1) produces 
formulae for Euler's constant. However, it should be noted that the following 
formula for Euler's constant, valid for all real x > 0, is actually equivalent to (1), 
namely, 



7r 1 7T log 1 2 sin(7rx)| x 



- 7 =-loga; + — r+ —Tin : >. 



r-2 



3 to 4vrx 2 2sinh 2 (vTx) ^fc(£; 2 + x 2 ) 



2k 7r log \ k A — x 



^ ( k 2 + x 2 )(e^ k - 1) - 2 2^ sinh 2 {irk) 

_ 2 .| e --^|^_,3|«). (17) 

To see the equivalence of (1) and (17), consider the partial fraction expansion of 
7r cot(-7rx) 



3 27TX ^ 



in the form 



, . . 1 1 7T COt(vTx) 



Ex 2 ^ 4/cx 2 

i =i k(k 2 + x 2 ) + (e^ k - l)(fc 4 -x 4 ) 

(Entry 8, Chapter 30 [4, p. 374]). Subtracting (18) from (1), we obtain (17). 



(18) 
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If we use the classical formulae [1, p. 259, 6.3.13, 6.3.17] 
9^(1 + ix) = — coth(7rx) — — and 9fa/>(l + ix) = —7 + ^ 



x 2 



^ k(k 2 + x 2 ) 

and substitute the latter in (17), after cancelling —7 from both sides there results 
a formula for ip(l + ix), valid for x > 0, analogous to formula (1) for ^(1 + x). 

The evaluation (17) was originally proved by Bruce Berndt, Jonathan Borwein, 
and Will Galway. They showed that the derivatives of both sides of (17) vanish, 
and that the right side of (17) tends to — 7 as x tends to infinity. Since (1) and (17) 
are equivalent, this argument provides a verification of Ramanujan's formula. The 
question remains as to how Ramanujan arrived at (1) in the first place. We may 
never know how Ramanujan arrived at some of his results, but in the case of (1), 
we are able to give a plausible answer to our question. 

Ramanujan was in the habit of expanding meromorphic functions into their series 
of partial fractions, and many such examples can be found in the unorganized pages 
of the second notebook. However, the following partial fraction expansion (Entry 
3 of Chapter 30 [4, p. 359]) has a particular relevance to (1). 

7T 2 CSC 2 (7TX) 7T 1 1 i 

e 2 ™ - 1 = 3x" " 2X 2 + 2^x3 " ^ (X + } 

^ 4kx ^ 27tx 3 

+ (e *rfc _ 1)(x 2 _ £2)2 - sinh 2 (7rfc)(x 4 _ £4) ■ (ig) 

Observe that (19) is simply a rearrangement of (9). One is tempted to speculate 
that Ramanujan, observing the presence of the derivative ip (x + 1), was motivated 
to integrate both sides of (19) to obtain a formula for ip(x + 1). Carrying out the 
anti-differentiation reveals a happy correspondence between the terms that arise 
and several of the terms appearing on the right of (1). To indicate just one ex- 

i i . 7T ^ log I k A — X 4 1 . . . 

ample, we see the mysterious — — : — 2 appearing as the anti-derivative 

k=i i'ftk) 

00 



of \^ . However, a few terms on the right of (1) remain unac- 

^sinh 2 (7rA:)(x 4 - k 4 ) 



counted for, of which the double series is the most glaringly obvious. We shall see 
that the double series arises fairly naturally in the course of proving (1). Accord- 
ingly, we now focus on completing this final task. 

Proof of Theorem. By the product rule for differentiation, we have 

d 7TCOt(7Tx) 7T 2 CSC 2 (7TX) 7T 2 COt(7Tx) 

dx e 2 ™ - 1 ~ e 2 ™ - 1 2sinh 2 (7rx) 

and 

d it log |2 sin(7rx)| 7r 2 cot(7rx) tt . . . .. d 



+ — log |2sin(7rx)|— csch (ttx) 



dx 2sinh 2 (7rx) 2sinh 2 (7rx) 2 dx 
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Hence, we can rewrite (19) in the form 

,/. 7r 1 1 d ir cot(irx) d it log |2 sin(7rx) 
ib (x + 1) = I I - - H — - - 

VK ' 3x 2x 2 2ttx 3 rfi e 2 ™ - 1 dx 2sinh 2 (7rx) 

E4kx 2irx 3 
^ (e 2 * k -l) ( X 2_ fc 2)2 " si nh 2 (7rA:)(x 4 -fe 4 ) 

— — log |2sin(7rx)|-^-csch 2 (7rx). 

But, 

2dx CSC " fa ( e 2nx _ if ~ \dx J e 2nx - 1 ~~ \dx J ^ 



X>2tt*)V 



-2-nkx 

— an rv j c 

k=l 



Using this in (20) and integrating both sides with respect to x, we obtain 
7r , 1 1 7rcot(-7rx) 7r log |2 sin(7rx)| 



3 2x 4vrx 2 e 2 ™ - 1 2sinh 2 (7ra;) 

E2/c 7T log | k A — x A I 

^ pTij^T^j " 2 ^ sinh 2 (7rA:) 

- / log|2sin(vrt;)|^(-27rA:) 2 e- 2 ^dt; + C, 

J* k=i 

where C is a constant of integration to be determined. We now let 

Zoo 
log|2sin(vrt;)| ^{-2^kf e' 2 ^ dv 
k=i 

= / log|2sin(7r(x+n))| V(-2vrA;) 2 e- 2 ^ (a;+u) dn. 

Jo 



k=i 



By means of the Laplace transforms 

r°° n f°° k 

I e~ kt sm(nt)dt = ^ an d / e~ kt cos(nt)dt 

Jo k + n J 



k 2 + n 

and the well-known Fourier series 

E° cos(27r/cy) ln . . .. 
^ = -log|2sin(7n/)|, y real, 
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we find from (22) that 

00 /»oo 

S(x) =(2tt) 2 J2k 2 e- 27rkx / e- 27Tku log\2sm{Tr(x + u))\du, 
k=i Jo 

= - (2vr) 2 J2 k 2 e~ 2 ^ r e~ 2 ^ J2 COS(27rnX + ^ du 

k = l Jo n=l n 

= - 27T f; k 2 e~ 2 ^ e~ kt jr COS(2?rnX + nt) dt 

k=l ^° n=l U 

=2tt ^ k 2 e~ 27Thx [°° e~ kt ^ sin ( 27ma sin( nt) - cos(27rnx) cos(ni) 

fc = l ^° n=l 



k=l n=l " 

^TVe" 2 ^ f V sin ( 2?mx ) [°° e ~ kt sm(nt) dt 



=2, ± e~ 2 - ( A; 2 £ - fc 3 £ ) , (23 ) 

^ \ ^ k 2 + n 2 ^n(k 2 + n 2 ) v 7 

fc=l \ n=l n=l v y / 

where the inversions in the order of summation and integration are justified by 
absolute convergence. Substituting (23) for the definition of S(x) in (21), we have 

. . 7T 1 1 7TCOt(7Tx) 7T log |2 sin(7Tx)| 

^(x + l)=C+-logx + --— ^+ I + 



3 6 2x 4vrx 2 e 2 ™-l 2sinh 2 (7rx) 

E° 2k 7r ^ log I k 4 — x 4 I 

^ ( e 2^-l)(fe2_ x 2 ) " 2 A. sin h 2 (^) 

- 2^ V e - 27rna; C A: 2 V sin ( 27rnx ^ - A: 3 V cos ( 2?rnx ) ^ 

S V £1 fc2 + n2 ~i "(* 2 + "'V " (24) 

By letting x — > +00 in (24), we can evaluate the integration constant, C. In fact, 
C must equal zero, as a comparison of (24) and (1) reveals. Since most terms in 
(24) vanish asu +00, it suffices to show that 

(co 1 1 4 ,41 \ 
3 2 smh (7rfc) J 

Set x = N + 1/2, where iV is a positive integer, and write 

log|x 4 — A? 4 | / V s 1 l°g |a; 4 — k 4 \ 



f-? sinh 2 (7rfe) If-'. sinh 2 (7rfe) 

fc=l v y \fc<V^ ^<k<x k>xj v 7 

E41ogx ^ log(l — k 4 /x 4 ) 
sinh 2 (7r/c) sinh 2 (7rfc) 
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+ (4x log x) O (V 27r ^) + O (V 3 ™/ 2 ) 

= (^>(g-^ + ^-"i) 

+ 0(x" 2 ) + O (e~*^ 

oo 1 

= (41ogl) g5S?w +0(l " 2) ' <26) 

Second, by Stirling's formula, 

V>(x + 1) = logx + 0(l/x), (x^+oo). (27) 
Employing (26) and (27), we see that 

log |x 4 — k 4 \ 



^{x + l)-\^ogx + ^Y^ 



2 ^ sinh {irk) 

(OO \ 

as x tends to +oo. However (Berndt [2, Prop. 2.26]), 

7T 1 7T / 1 1 \ , 

1 - - + 2tt > ^ = 1 - - + 2tt =0. 29 

3 ^sinh 2 (vrA:) 3 V6 2vr y V ' 

Using (29) in (28), we find that 

^ + 1) _| logI+ |g^J. 0(1/l) , (30) 

as x tends to +oo. Thus, the limit (25) holds, and the proof of the theorem (1) is 
complete. □ 
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